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Pages 34 and 35: The equations for and the coefficients 

I * » * ^N. should be modified 
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METHODS OF ANALYZING WIND-TUNNEL DATA 
FOE DYNaKIC flight CONDITIONS 
By C* J* Donlan and I. G. Recant 



SUMMARY 



The effects of power on the stability and the control 
characteristics of an airplane are discussed and methods 
of analysis are given for evaluating certain dynamic char- 
acteristics of the airplane that are not directly discern- 
ihle from wind-tunnel tests alone. Data are presented to 
show how the characteristics of a model tested in a wind 
tunnel are affected hy power. ' 

The response of an airplane to a rolling and a yawing 
disturbance is discussed, particularly in regard to changes 
in wing dihedral and fin area. Solutions of the lateral 
equations of motion are given in a form suitable for direct 
computations. An approximate formula is developed that 
permits the rapid estimation of th^ accelerations produced 
during pull-up maneuvers involving abrupt elevator deflec- 
tions. 

INTRODUCTION 



Some time ago, the NACA undertook an investigation to 
determine the flying qualities of a low-wing, pursuit mon- 
oplane. • This airplane (fig. 1), like many of its type, 
was found to possess several undesirable flying character- 
istics. Accordingly, a number of modifications to the 
airplane were recommended. In order to study the effects 
of the proposed modifications and in order to investigate 
the. reliability of wind-tunnel: data for estimating the 
behavior of the full-scale airplane, a scale model' of the 
airplane, with and without the proposed modifications, 
was- tested in a wind tunnel. '■ Tests were made with and 
without the propeller operating in order to study the ef- 
fect of power on the characteristics of the model. 



2 



KACA Technical Note No. 82-8 



During the course of these w'ind— tunnel tests, several 
interesting results were ohtained concerning the effects of 
power on the characteristics exhibited hy the model. In 
the evaluation of the data, certain methods of analysis 
were employed that were found useful in estimating the 
dynamic "behavior of the airplane, particularly as regards 
the effects of the proposed modifications. These subjects 
are treated in the present paper. It is fielt that this 
, paper will serve as an aid to other wind-tunnel investiga- 
tors seeking to evaluate the flight characteristics of an 
airplane from wind-tunnel results. 

The first portion of the paper is concerned with the 
effect of power on certain of the wind-tunnel character- 
istics exhibited by a model. &raphs are presented that 
emphasize the importance of obtaining wind-tunnel data 
with powered models. The second portion of the paper dis- 
cusses the lateral motions of the airplane. Concrete ap- 
plications are included showing how the solutions to the 
lateral equations of motion were used to estimate the 
change in dynami.c flight characteristics likely to result 
from the use of the recommended modifications to the air- 
plane. Finally, the longitudinal motion of the airplane 
is considered. An approximate formula is developed that 
was found useful in estimating the accelerations and the 
stick forces . likely to develop in an airplane during an 
abrupt pull-up maneuver. 



THE IMPOBTANCB OP POWEEED MODELS IN WIND-TUNNEL TESTS 



Although the influence of the slipstream on the char- 
acteristics of a:n airplane has been appreciated for many 
yearst wind-^tuhnel tests with propellers operating have 
not. been the usual procedure.- The frequency with which 
undesirable flying characteristics appear in modern air- 
planes has led- to some apprehension concerning the relia- 
bility of wind-tunnel data obtained without powered models. 
The particular testing technique employed in the present 
invest igat ion will be reserved for detailed discussion in 
a future paper. It will suffice to say here that the 
■power-on test- procedure should be such that data' may be 
obtained for thrust conditions of the model corresponding 
to those of the full-seale airplane at the lift coefficient 
and the flight attitude under consideration. 
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Figure 2 is a plot o f the* .el evato r angle required for 
trim against the indicated air speed at whicli trim occurs. 
These curves were evaluat ed f rom data secured from (l) 
wind-tunnel tests.witht no propeller on the model, (2) ac- 
tual flight tests with the full-scale airplane, and (3) 
wind-tunnel tests with a powered model. The shapes of 
the curves are indices of the static ' sta*bility (see later 
discussion of longitudinal motion) displayed by the air- 
plane in flight and "by the models in the wind tunnel. The 
agreement hetween the flight results and -t-he results ob- 
tained with a powered model is good:,part icularly for the 
flap-up condition. The larger discrepancy existing for 
th-e flap-down condition is partly attributable to the 
rather large tab deflection (11° nose up) used on the 
full-scale airplane" for this test.- The related variation 
of the elevator ' stick force required for trim with indi- 
cated air speed is not presented here because the model 
data and the flight data were not direct ly compairable . . 
owing to the absence of a trim tab on the' "model. • "In the 
evaluation of the information obtained from wind-tunnel 
tests, certain methods of presenting the data have been 
found to facilitate the analysis. The' result ing' pitching 
moment acting on an airplane trimmed for. steady equilib- 
rium flight is zero. This condition is true whether or 
not the airp-lane is flying 'yawed. As far as the pilot is 
Gbncerned, the dihedral effect of the wing manifests it- 
self as a pure rolling moment about an axis lying in the 
plane of symmetry of the airplane. In wind-tunnel in-' 
vestigat ions , oh the other hand, the rolling moment is 
usually measured about an axis coincideht with the tunnel 
axis. Consequently, when the model is yawed relative to 
the tunnel axis,- the measured rolling moment about this 
axis will contain a component of th.e pitc^iing moment act- 
ing: on the model unless the model is trimmed for zero 
pitching moment. In thie comparison of the slopes of 
different rolling-moment curves obtained* in yaw tests, it 
is therefore necessary to know whether the measurements 
wiere made with the model trimmed in pitch. - 

If the system of axes discussed in the subsequent 
section on lateral motion is employed for the presentation 
of wind-tunnel ' data , the necessity* -for cont inually check- 
ing the trim of the model during yaw tests is avoided. 
In addition, the data in this form can be used directly 
in the evalxiat ion of the stabll;ity del'iv-at ives necessary 
for- the solutions of the equations of motion. * 
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In figure 3 are plotted for various power conditions 
the rates of change with angle of yaw of :the yawing-^ moment 
coefficient dCn/i^» the : rolling-moment coefficient 
dC-^/d\|;v and the lateral- force coefficient dCy/d^t^. • The 
decrease in dihedral effect dC^/d\l/ with the application 

of power should be ohserved. In the condition presented 
(flap up), the decrease , in dihedral effect with the applica- 
tion of power is not great, but, in the flap*down condition 
and at high lift co ef f icient s.i the effect of power may be 
critical • 

An inspection of figure 4 will suffice to illustrate 
the effect of the slipstream on the. charact erist ics of a 
control -surf ace placed within its boundaries. The nega- 
tive yawing moment that exists with neu-tral rudder when 
the propeller, is operating should be observed. The moment 
is apparently caused by the rotation of the slipstream. 
It is, however, recognized that the slipstream rotational 
characteristics exhibited by the model in this test are 
probably nat identical with the slipstream rotational 
characteristics that would be obtained in flight with the 
full-scale airplane, even under similar thrust conditions. 
The results do indicate, nevertheless , that with the pres- 
ent trend toward greater and great er -power , the effect of 
this slipstream rotation may be a critical factor in tail 
design. 

Figures 2, 3, and 4 have been included to illustrate 
the effects that power may have on the stability . character- 
istics of an airplane and on the effectiveness of the hor- 
izontal and. the vert ical tail surfaces. These effects are 
difficult to . estimate and consequent ly recourse must be 
had to ' wind-tunnel tests of a model equipped with running 
propell€|rs« Wind-tunnel tests of a powered model require 
more care and are more expensive to perf o.rm .than tests of 
the conventional type of model. By a judicious choice of 
tests, however, it is possible. to secure the information 
of greatest value with a minimum of time and expenditure. 

..THE-LATESiiL MOTIONS OP. T5E AIRPLANE. 



Certain lat eral-r stability and control characteristics 
of an airplane can be directly determined from wind-tunnel 
tests of a model. These characteristics usually include 
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(1) the directional stability (fig. 3(a)), (2) the dihe- 
dral effect (fig. 3(h)), and (3) the static effectiveness 
of the control surfaces (fig. 4). This information is in- 
dispensable to the designer if he is to proportion his • 
airplane properly and is also necessary for the evaluation 
of the stability derivatives from which the dynamic be- 
havior.of the airplane may be predicted. 

The dynamic flight characteristics that an airplane 
will display are not, however, readily obtained. The be- 
havior of an airplane subjected to a lateral disturbance 
is complicated by the coupling of the ensuing rolling and 
yawing motions. The coupling of these two motions makes 
it impossible to estimate the effectiveness of the ailer- 
ons and the rudder in flight from static tests of the con- 
trols, alone. A sound evaluation of these factors can be 
made. only if a knowledge of the motions produced by the 
controls is available. 

The angle of bank produced, for example, in 1 second 
by the total deflection of the ailerons might be taken as 
a measure of the effectiveness of the ailerons in produc- 
ing a sudden roll. If, on the other hand, the'airplane 
is rolled slowly, a greater yawing motion being permitted 
to develop, the result may be to produce a rolling motion 
opposed to that generated by the ailerons and of such 
magnitude that the effect of th^e ailerons is completely 
nullified.. The airplane may even ro'll against the ailer- 
ons. 

In actual flight the pilot might counteract this ad- 
verse rolling tendency by coordinating the rudder with the 
ailerons. In the present analysis, however only the in- 
herent dynamic characteristics of the airplane are con- 
sidered. It would be of interest, then, to know the varia- 
tion of the angle of bank with time after small aileron 
deflections are applied. . For similar reasons, the effec- 
tiveness and the sensitivity of the rudder can be most 
carefully judged if the angle of bank and the angle of yaw 
produced in a definite time interval by a small rudder de- 
flection are known. Information of this type is particu- 
larly valuable in e st imat ing t he effect of various combina- 
tions of vertical fin area iand wing dihedral on the dynam- 
ic behavior of an airplane. 
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Apart from actual flight tests, pro^oably the hest 
method of estimating the dynamic qualities of an airplane 
is by evaluating the dynamic equations of motion, the data 
obtained from w.ind-tunnel tests of a powered model being 
used. 

Assumptions- and symbols,*- The assumptions generally 
made in the study of airplane stability are made here. 
The most important of these assumptions are: 

1. The air forces* and the moments resulting from 
displacements of the airplane relative to its steady con- 
dition of flight are proportional to the displacements or 
to their rates of change. 

2» The components of moment due to different compon- 
ents of the motion are directly additive. (For example, 
the rolling moment due to combined* rolling and sideslip-, 
ping may be computed as though the rolling and the side- 
slipping had occurred separately.) 

The axes used in specifying the moments, the angular 
velocities, and so forth are fixed in the airplane and 
move relative to the earth and the air. The X axis, 
passing through the center of gravity of the airplane, is 
in the plane of symmetry and is so oriented that it points 
into the relative wind when the airplane is flying stead- 
ily in unyawed flight. Also, the aixes form a conventional 
orthogonal system intersecting at the center of gravity. 
The Z axis points directly downward in the plane of sym- 
metry and" the Y axis points aldng the direction of the 
right wing. The motions discussed iare those, of the moving 
axes relative to the undisturbed air with the exception of 
the angle of bani, which is measured' f:i*om the horizontal. 

• The symbols used in the following analysis are de- 
fined in the appendix. • * 

Although the axes change their orientation in the 
airplane with different lift coefficients and probably,: 
never coincide with the axes of the principal moments of 
inertia, the corrections in unstalled 'flight are small 
and have been neglected, as have the products of inertia. 

Equations of motion,- If the airplane is considered 
capable of motion in all degrees of lateral freedom, the 
equations of motion with deflected controls (neglecting 
the small side forces developed by the deflected controls 
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and by the rolling and the yawing. velocit ies) may be writ 
ten: 



Also 



dp 
dt 



- = pLp + rLy + pLp + 6Lg (in rolling) 



dr 
dt 



^ = pNp + rNy + pNp + 6Ng (in yawing) 



pping) 



dt ^' dt • ^ Uo 



(1). 



In order to solve for any of the variables, it is 

necessary to integrate this system of linear simultaneous 
equations* For reasons that will be apparent later, it 
is convenient to solve the system. of equations (l) for 
separate unit magnitudes of the control disturbance terms 
6L5 and 6N5; that is, one set of solutions is obtained 
by letting 6Lg = 1 and 6Ng = 0, and another set of 
solutions is obtained by letting 6Lg = 0 and. 6Ng = 1. 

The unit disturbances are assumed to be instantly applied 
at zero time and to remain constant thereafter. In order 
to distinguish the separate solutions, the subscript L 
will be applied to the solutions obtained when 6Lg = 1 

and 6N5 = 0; whereas the subscript N will be applied 

to solutions obtained by letting 6Lg = 0 and 6Ng = 1. 

Thus, pj^ represents the rolling velocity resulting from 

the application of a pure rolling disturbance of unit 
magnitude. 

If the symbol D is substituted for d/dt and if 
&Lg = 1 and 6Ng = 0, equations (l) may be rewritten in 
the. following form: 
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p (D-Lp) - rLj. - pLp = 1 
p (-Np) + r(D-Nr) - fiNp = 0 

p L + Dr + D ("d - -^^ p 



= 0 



(2) 



It can "be shown that the solutions of equation (2) 
are of the form: 



X, t 



Aot 



X,t 



PL = + Pl,^"'" + PLs^ ^ Pl^^ ^ + PL4«, 



^4t 



(3) 



where 



p; resultant rolling velocity due to unit 

rolling disturhance, that is, '5L§ = 1 



Pl "Pt 
^0 -^4 



constants that depend only on values of 
stability derivatives and on type of . 
disturbance involved 



••■^4 J^oots of stability equation F(D) = 0 



where 



F(D) = 



D-L, 



-s. 



-L, 



D-K, 



-L 



• The constants in equation (3), together with the con- 
stants appearing in the expressions f.or all the remaining 
components of motion (including the solutions for the unit 
yawing disturbance 61T§ = 1)» have been evaluated in terms 
of the stability derivatives and the roots of the stability 
equation. The solutions of equation (2) are tabulated in 
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the appendix- in a form suitable .fo r computation. The 
jstability equation is also discussed and an alternative 
form of equation (3) is- given fo.r use when the solution 
of P(D) = 0 includes conjugate complex roots of the form 
a ± ih. 

The expressions presente.d in the appendix were eval- 
uated hy applying the operational mathematics of Heaviside 
to equation (2). For the theory of operational methods 
the reader is referred to a standard text on the subject, 
for example, reference 1. Specific applications of the 
Heaviside treatment to other problems in airplane dynam- 
ics may be found in reference 2. 

After the complete unit solutions have been obtained, 
time histories of the motion caused by the unit disturb- 
ances can be plotted with very little additional calcula- 
tion.. Because of t he. 1 inear ity of the equations of motion, 
the unit solutions may be compounded in any arbitrary man- 
ner. If, for example, 6^15 represent's the rolling ac- 

a 

celeration created by the applied aileron rolling moment 

and S^No represents the accompanying yawing accelera- 

^ a 

tion, then, at time t, 

^t -^aH^ (f^X + Sa^S.^ - (4). 

. where 

resultant angle- of bank after t seconds 

(^Ljt ^'^Sle. of bank after t seconds due to a 
^ unit rolling disturbance 

(^Sj angle of bank, after t seconds due to a 
^ unit yawing disturbance. 

As. in the case of the unit disturbance terms, the 
actual disturbances, 6«Lc and 6«Nc , are assumed to 

be suddenly applied at zero time and to remain constant 
thereafter. 
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A ppl icat ions*- The practicability of using the solu- 
tions to the lateral equations of' motion. is best illus- 
trated by concrete applications. 

The airplane under consideration, as mentioned in the 
introduction, exhibited certain undesirable flying charac- 
teristics in flight tests: The dihedral effect was un- 
desirably low at slow speeds- with the flaps down (fig. 3(b) 
indicates satisfactory dihedral characteristics with the 
flaps up over the range considered), and.it was impossible 
to raise a wing by use of the rudder alone. In an' effort 
to improve the lateral flying qualities of the airplane, 
it was proposed to increase the wing dihedral. Because 
the lateral characteristics of an airplane depend not only 
on the absolute amount of dihedral but also, on the relative 
amount of weathercock stability present, it was considered 
necessary to increase the vertical tail area as well as 
the wing dihedral. The ailerons were unchanged, but the 
rudder was so modified as to improve its hinge-moment 
characteristics. These modificat ions are shown on figure 
1. 

Wind-tunnel tests' were made with a model equipped for 
power-on tests with and without the proposed modifications. 
The data from the comparative tests indicated that con- 
siderable improvement should result from the incorporation 
of the modifications on the full-scale airplane. For all 
flap and power conditions the dihedral effect dC^/d\!/ re- 
mained positive, the index of weathercock stability 6.C^/d^ 

remained negative, and the static characteristics of the 
rudder appeared satisfactory. Wind-tunnel tests of a 
model, however, provide no direct information pertaining 
to the dynamic flight characteristics of the full-scale 
airplane. Instances have occurred in which the incorpora- 
tion, of similar modifications on a full-scale airplane has 
affected the control characteristics of the airplane in an 
adverse manner, particularly* at high speeds, in spite of 
the favorable static characteristics indicated from wind- 
tunnel tests. Accordingly, it was decided to investigate 
the response of the airplane, with and without the modifi- 
cations to the aileron and the rudder controls, by evaluat- 
ing the solutions to the equations of motion.- 

411§.?l?.*^I_?.o^E^rQLl'~ particular problem considered 

here, undesirable aileron control characteristics are like- 
ly to be manifested in the form of aileron "heaviness" or 
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stiffness at high speeds. The suhsequent analysis, how- 
ever, is perfectly general and is in no way limited to 
this high-speed condition. This particular condition is 
treated, laerely as an example, to indicate the general 
method of procedure. Because the ailerons are identical 
on hoth the original and the modified airplanes, aileron 
heaviness can be physically interpreted as an increase in 
stick force resulting from the increased aileron deflec-- 
tion necessary to reproduce a given rolling maneuver with 
the modified airplane-. This . int erpretat ion of aileron 
heaviness suggested the fo llowing method of analysis: 

(a) On the assumption that the original airplane 

(airplane A) was flying in steady high-speed flight, the 
angle of hank generated in 5 seconds hy a small aileron 
deflection 6^ computed. 

(h) Then, the aileron deflection 6^ + necessary 
to reproduce the identical maneuver with the modified air- 
plane (airplane B) was calculated.. 

The time for the maneuver (5 sec in this case) is 
somewhat arbitrary. It should be of sufficient duration, 
however, to permit full development of the secondary roll- 
ing effects introduced by the induced yawing motion. 

The magnitude of the increment A5a is a direct 
measure of the additional stick force required to perform 
the maneuver (because the ailerons are identical) and may 
be used as an index of aileron heaviness.. If A5g^ is 

large and positive, the aileron stick forces on airplane 
B may be too large to be acceptable and modification of 
the aileron itself may prove desirable. 

The stability derivatives of airplanes A and B 
were evaluated from. the power-on wind-tunnel tests of the 
two models and from data in references 3' and '4.. These de- 
rivatives are tabulated in table I, together with other 
information necessary to evaluate the solutions to the 
equations of motion. 

In accordance with equation (4). the angle of bank 
assumed by airplane A in t seconds after the applica- 
tion of an aileron deflection 5-.;. is 
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In order to evaluate the unit solutions and 

it is necessary to oibtain the roots of the stability equa- 
tion F(D) = 0. If the appropriate values from table I 
are substituted into the expression for F(D) given in the 
appendix, 

F(D) = D^ + 20.4555D^ + 52.78841?^ + 347.8242D + 5,4S7 60 = 0 

With the use of the procedure outlined in the appendix 
for solving quart i c equat ions , the roots of this equation 
were determined to be: 



\ 


= -0.01567 




K 


= -18.6230 






= -0.908424 + 


4.2199i 




= -0 . 908424 - 


4. 21 99 i 



Substitution of the appropriate roots and derivatives 
in t*h:e expression for the angle of bank , " give'h in the 

•appendix, yields 

i . — 0.01 56 7 t -18,623t 

0j = 3.3471 3.3497e + 0.002876.e 



+ «0 •.167 7 8 _0 i.1 p 3 1^ 

V 1741.79 - 2270.451 / 

+ /^-0 > 46778 >- 0>16310 93i'\ 

V 1741.79 + 2270. 45i / 



(-0.9034 + 4,2199i)t 



e 



(-0.9084. - 4,2l99i)t 

e 



If the last two terms are combined in accordance with 
the transformation formula in the appendix, the final ex- 
pression for ^® written: 

, -0.015 67 1 — lS.623t-;-. 

= 3.3471 - 3.3497e . + 0.002876e 

- 0.00034623e'^^'^^^^^ cos 4.2199 (t + 0.1354) (5) 
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S imilarly , 



i — 0.01567t — 18«623t; 

0^ = 11.5982 - 11.62275e + 0.00038854e 



+ 0.043548e ^'^^^^^ ^os 4.2199 (t + 0.23305) (6) 



Practical solutions are most conveniently obtained 
"by graphical addition and subtraction of the component 
parts of the motion. It can be seen that a subtraction 
involves the small difference of relatively large quanti- 
ties. For this reason it is necessary to retain as many 
absolute figures as possible in the evaluation of the 
roots of the stability equation and in the evaluation of 
the individual components of the motion. 

Equations (5) and (6) are plotted in figure 5 (a) 
for airplane A. The corresponding expressions for air- 
plane B are plotted in figure 5 (b)." 

Prom figure 5 (a) the unit solutions for and 
02J after 5 seconds are 

0.25 radian 
0.86 radian 



Por a 1° aileron deflection from table I, 
6aL;> = 1.54 -oer second. 

and 

SalTfi^ = 0 

Hence, for airplane A, 

{0) 5 = (1.54) (0.25) = 0.385 radian = 22.1^ 

The to tal' a ileron deflection necessary to bank air-* 
plane B 22.1^ in 5 seconds is given by 



' T 

^ 5 sec 



5 sec 
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Sa = 



where 



^6 



mk. ^ 



Z 

For airplane Lg =1.54 per degree, Ng = 0, 



a 



and = 0.26 radian; hence. 



0.385 
(1.54),. (0.2 6) 



6a = ."^r , . = O.S6^ 



or 



= -0.04° 

In view of the fact that the ailerons on the two air- 
planes possess the same hinge-moment characteristics, it 
can "be concluded that the aileron^- stick- fo rces developed 
during maneuvers will at least "be no greater (or heavier) 
on the modified airplane than on the original airplane, 

^^i^der _c_o^nt^ro 1.-- The sensitivity of the rudder con- 
trol, particularly at high speeds, may be judged by the 
angle of bank and the angle of yaw generated, after a def- 
inite time interval, by the yawing moment impressed by 
the rudder. Here again the time must be of sufficient 
duration to permit the interaction of secondary rolling 
and yawi-hg effects. 

Airplane A w^s assumed to be flying level at. high 
speed and to be suddenly, subjected to a yawing moment 
impressed by a given rudder deflect ion. The angle of bank 
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of and the angle of; yaw \l/ general ed after 5 seconds, 
were calculated. Under similar circumstances, the angle 
of bank 0 + and the angle of yaw \1/ + A\i/ developed 

\>y airplane B were similarly calculated. The magnitudes 
of and A\i/ were taken as a measure of the relative 

sensitivity of airplanes A- and B to their respective 
rudder controls. The computations follow. 

From figure 5 and table I"-, 

Ai rpla ne _A A lrip lane B 
radian .86 .98 



radian - 0.25 0.26 

5 sec 



5 sec 



6^L;. , per sec^ ,308 .410 

r 



8j,Ng^, per sec^ .549 . 652 



If a form similar to that of equation (4) is used, 

and the angle of bank of airplane A, 5 seconds after the 
sudden application of a 1^ rudder deflection, is 

= (-0.308) (0.25) + (0.549) (O.86) 

o sue 

= 0.395 radian = 22.6^ 
For airplane B, 

^ ^^r. ^ (-0 .410) (0.23) - (0. 652) (0.98) 

= 0.532 radian =30.5^ 

and 

. A0 = 8^ ~ 
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The angles of yav? generated "by the unit disturbances 
were obtained in a •similar manner: thus, 



In the maneuver investigated, the unit solutions after 
5 seconds were as follows: 



Airplane A Airplane B 

radian ------- 0.04 0.04 

5 sec 

radian - -- -- -- ,30 .26" 



^5 sec 
For. airplane A, 

sec ^ (-0.308)(0.04) + (O . 549) (O . 30 ) 

= 0.153 radian = 8.8^ 

and, for airplane B, 

(\|/) = (-0.410) (0.04) + (0. 652) (0.26) 

5 se c , 

= 0.154 radian = 8.8° 



0 r 



It appears that, for equal rudder deflections, the 
modified airplane will tend to generate more bank than 
the original airplane. The increase in the wing dihedral 
and the modified fin and rudder may restilt, then, in a 
rudder control that will be slightly more sensitive at 
high speeds than the rudder control on the original air- 
plane. 

Einge-moment measurements are, of course, necessary 
to determine whether the resultant rudder-pedal forces on 
the modified airplane will be greater or less than those 
on the original airplane. .The hinge-moment characteristics 
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Por airplane A, R = -0.297; for airplane B, E 
= -0.150. The resultant pedal forces, on the other hand, 
depend not only on the hinge-moment .characteristics of 
the rudder but also oh its floating characteristics and 
on the particular deflection required to perform a stii:- 
ulated maneuver. The resultant hinge moment per degree 
of rudder deflection may be expressed as follows: 



The quantity (^^q q theoretically a con- 

r 

stant for any given tail arrangement but actually it is 
very critical to interference effects at the tail and 
fluctuates considerably. For small angles of yaw (±5°), 

\ 'S\ir/C 0 practically zero for both airplanes A 

r 

and B; and the quantity R may, therefore, be taken as 
a measure of the pedal forces for a given rudder deflec- 
tion. 



THE LONGITUDINAL MOTION OP THE AIRPLANE 



When the longitudinal, stability of . an. airplane is 
discussed, the characteristic usually referred to is the 
"static" longitudinal stability. If static longitudinal, 
stability exists, the dynamic stability characteristics 
are of minor importance (reference 5). 
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The usual index for static longitudinal stability is 
the rate of change of pit ching-raoment coefficient with 
lift coefficient dC^^/dCL or some quantity proportional 
to it. In flight the most convenient method of evaluating 
the amount of static stability present is to measure the 
elevator angle required to trim the airplane at various 
speeds; the slope, d6^/dV or dg^/da, is an index of 

the degree of static stability possessed by the airplane. 
The significance of the ratio d^e/^^* methods of eval- 
uating it for power-off and windmilling conditions, and 
suggested design values are discussed in reference 6. The 
necessity for power-on wind-tunnel tests for securing., the 
effect of power has already been discussed (fig* 2). 

If tunnel data are available from which the floating 
angle of the elevator may be caluclated for any lift coef- 
ficient, the elevator stick force required for trim at any 
lift coefficient can.be calculated from- the following 
f 0 rmula: 



/ 



\ dx 



^3 C 



\d6 , 



(7) 



where 



A5 



P 

X 

o 



stick force for trim, pounds 

linear travel of top of control column^ feet 

difference, in degrees,, between elevator 
. angle. required -for trim. at lift coef- 
ficient under consideration and free- 
floating angle of elevator 



The term A6 ^ may fluctuate considerably with power. 



The slope of the curve relating the variation in stick 
force with forward speed dp/dV depends on the speed at 
which the airplane is tritamed for zero stick force and, 
consequently, depends on the initial setting of the trim- 
ming tab. In the comparison of the change in the slope of 
the stick-force curve resulting from modifications to the 
elevator, care should therefore be taken to orient the 
elevator t rimming-^ tab settings so that zero stick force 
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always occ-ars at the same speed. Otherwise, a superficial 
examination of the curves of stick-force variation with 
air speed may lead to incorrect conclusions concerning the 
effectiveness of the tail surfaces or the stability of the 
airplane^ 

In addition to providing a means of trimming the air^ 
plane in steady flight, the elevato-r must he capable of 
changing the airplane flight path. The rate at which this 
change is accomplished in a quick pull-up can, in a sense, 
"be interpreted as a measure of the effectiveness of the 
elevator in maneuvers. 

It is convenient to take the rate of change of the 
maximum normal acceleration per unit of elevator deflec- 
tion as an index of elevator effectiveness in maneuvers. 
Care must "be exercised, however, in interpreting this 
index. Although it is essential in a pursuit airplane to 
design an elevator sufficiently powerful to maneuver the 
airplane to the maximum lift coefficient of the wing, it 
has been found very undesirable if this condition is ful- 
filled with a minimum amount of elevator deflection. As 
discussed . in., reference 6, satisfactory static stability 
characteristics require the quantity dfi/da to have a 
value around 0.5. In airplanes that required consider- 
ably less stick travel to trim the airplane over the 
angle-of-at tack range inadvertent stalling has frequently 
occurred in accelerated maneuvers. The optimum value for 
the rate of change of normal acceleration per unit of 
elevator deflection is. therefore conditioned by the re- 
quirements of satisfactory static longitudinal stability. 

In the fo llowing .sect ion a simplified formula is 
developed . that permit s. the . rap id e st imat ion of the normal 
accelerations developed in abrupt pull-up maneuvers. 

De velopm ent of a si mp lified formula^for normal ac- 
cele ra tion s prod\iced in abrupt p ull-TOs To a first 
order of approximation, the equations of. motion in the 
plane of symmetry involving a disturbance in pitch may be 
written as follows: 

~ = ^-"hi + + QXq - ge (a)^ 



dt 

dw 
dt 



= uZ^ + wZw + q (Uo + Zq)- gSoB + ^e^^ 



y (8) 



dt 



uM^ + wM^ + qMg + 6eMs (c) 

e 
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Also 




where 9^ is the initial angle the X axis makes with 
the horizontal :axis*" 

The axes in which equations (8) are expressed are 
similar to- those defined in the section on lateral motion 
except that the X axis is inclined at an angle 9 to the 
relative wind. 

In general, this system of equations must "be solved 
hy methods analogo-cis to those used in obtaining the solu- 
tions to the lateral equations of motion. The normal ac^ 
. • celeration is given by* the expression dw/dt - qUQ • 

These equations have been solved for the airplane, 
the characteristics of which (evaluated, so far as pos- 
sible, f rom wind- tunnel tests on a powered model) are 
presented in table II. The curve of the normal accelera- 
tion resulting from a 1^ upward elevator deflection is 
plotted in figure 5 together with other components of the 
motion* 

The variation of the components of motion shown in 
figure 6 is typical of an abrupt pull- up maneuver at high 
speed. The following facts are apparent from this figure: 

1, During the time the airplane takes to attain 
maximum acceleration, the velocity V remains sensibly 
constant and is equal to Uq i the equilibrium velocity. 

Accordingly, u = du = 0» and equation (8a) can be ne- 
glected. 

2, In the vicinity of the maximum acceleration, the 
.time rate of change of pitching velocity dq/dt is ap-* 

proximately zero, that is, q is approximately constant, 

Z. In the vicinity of the maximum acceleration, the 
acceleration component dw/dt is almost zero and the ac- 
celeration thereafter is given almost entirely by the 
product qUQ. 
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If Z§^, Zq, and §6^ (all of which are small) are 
neglected and the substitutions clZq^ = wZ^ and wM^ 
= aU(x are made, equation (8), for an abrupt pull-up ma- 
neuver, can be simplified to the following form: 

qUo + aZa =0 ^ 
qMq + alia = -Mq 



where 



6 



The solution for qU^ is 



where 



°- \da/ 2 m 
"a - \ia) 2 



pSV c 



In the expression for Mq , the slipstream factor P 
normally has a vaitie "between 1 and 1.25 for high-speed 
flight. The tail efficiency factor t]^ is always less 

than 1 and is generally ahoufO.S. The value of the prod- 
uct ■ Yt\^ is therefore always about unity. 
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If the airplane is assuined to "be flying level before 

the pull-up maneuver, mg =. Cj^ipSV^ = C^ipSUQ^. * If the 

2 ' 



suhst itut ibns = V, jjl = r , and ?rit - 1 are made, 

|PSI 

equation (9) can be reduced to the follov^ing simplified 
form: 



^ ^ gji I ^l^J^^^e 

fie -L r I fiCiV IL 

\ dCr c \da / S 



(10) 



where QUq/oq represents the change in normal accelera- 
tion per unit of elevator deflection.;. 

Formula (10) gives a value of the normal acceleration 
produced during an abrupt pull-up maneuver that closely 
approximates the raaximum acceleration which would be ob- 
tained by solving the more cJ^'mber some equations of motion. 
It is the acceleration that will be obtained if the ele- 
vator is instantly deflected to its final position and 
held in that position until the masiimum acceleration is 
reached. The formula gives le.ss accurat e ' resialt s when 

dCiu/dCL is small, that is, . - — <0.01. Several cal- 

ctxlations made for smaller values of. ^CrTi/dC^. gave re- 
sults, however, that were in error by less than 7 percent. 
Greater accuracy is to be expected at small values of Cj^ 

on accoimt of the "approxima.t ions made involving G. At 
high values of Cj^, these approximations introduce greater 

errors. 

ManeuverabiM The manner in which 

the normal acceleration produced iri a.- pull-up maneuver is 

affected by dCjjj/dCT and u is sho^vn in figure 7, An 

examination .of this f.igure ■ reveal s that, for airplanes 
with high wing loadings (large, y.) and low static* stability 
(small dCni/^iCL), the normal acceleration produced per 

degree of elevator is affected by a small change in the 
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static stalDility. In accelerated maneuvers this effect, 
as far as the pilot is concerned, will manifest itself in 
the form of increased stick forces if it is assumed that 
the pilot wants to produce a given acceleration regard- 
less of the degree of static stahility present in the air- 
plane. 

The index of static stability dCjjj/dCL for an air- 
plane with the characteristics given in taDle II is as- 
sumed to he -0.022. The curve of the normal acceleration 
produced in a pull-up from level flight at 448 feet per 
second is given in figure 6. The maximum change in the 
normal acceleration attained with an upward elevator de- 
flection of 1^ is about 84.5 feet. per second. 

Prom formula (lO) 

■ iHo ^ 3 4.4 -0.021 : ^ 

. c 0 ^ 0.10. V(34;4) (-0.022) + (-2.01)/ ^ " ^ ^ 

^0 



er 



degree and qU^ = 83.7 feet per second^, an error of less 
than. 1 percent. 

'If the static stability of the airplane is improved 
by moving the center of gravity forward 0.078c, the, ex- 
pression for the normal acceleration becomes 

• \ 0 ^:0.10 \(34.4) (-0.10) + (-2.07)/. .® 



If it is desired to load the airplane to its former 
additional load factor of 2.6g, the elevaitor movement re- 
quired with the higher ^Oj^/&Oj, would be 

e -1.37g 

Accordingly, the pilot would have to exert nearly 
twice as much sticK force to execute the identical pull- 
up maneuver because of the increased static stability of 
the aiirlane. The desirability of modifying the elevator 
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in order to change it s hinge-moment.-charact eristics thus 
partly, depends on the relationship, between the, static 
stability dedired and the magnitude of ..the st.ick .forces, 
acceptable in accelerated raaneuyers . 



CONCLUSIONS 

..This paper is intended to illustrate primarily how 
power affect s the charact efist ics o f a: inodel tested in the 
wind tunnel and ho\7 wind~!t urinel data may be used to esti- 
mate flying qualities not directly discernible from wind- 
tunnel tests. The analyses presented in this paper permit 
the following conclusions to be drawn concerning the meth- 
ods employed: 

(1) In the prediction of the flight characteristics 
of an airplane operating with power, considerable error 
may be introduced if wind-tunnel data from fests of a 
model not equipped with an operating propeller ate used. 

(2) In the analyses of wind-tunnel rolling-moment 
data, care should be taken, in the determination of the* 
dihedral effect, to allow for the contribution of the un- 

■ balanced pitching moment to the slope of the rolling- 
moment curve 

(3) The evaluation of the equations of motion permit 
estimates to be made of the relative effectiveness of the 
aileronst and of the' rudder- controls, particularly when 
changes in wing dihedral are involved as these effects are 
not readily discernible from static tests of the controls 
alone. The methods used have been found reasonably ac- 
curate, when wind-tunn-el data are available, and are not 
difficult to employ. 

.(4). The rate of change' of normal acceleration per 
unit of elevator defle.ction affords a convenient correla- 
tion between maneuverability, stability, and elevator- 
stick forces. The approximate formula developed for cal- 
culating the normal accelerations produced during an 
abrupt pull-up maneuver is simple to evaluate and has been 
found to yield reasonably accurate results, particularly 
at high speeds. 

langley Memorisil Aeronautical Laboratory; - 

National 'Advisory Committee for Aeronautics, 
Langley Field. Va. , July 22, 1941. 
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APPENDIX 
List of Symbols 

Uq velocity along axis in steady flight 

V velocity along flight path * 
n Velocity along X axis 

V s idesl ipping component of veio c ity 
w velocity along Z axis 

p angtilar velocity in roll 

q angular velocity in pitch 

r angular velocity in yaw 

0 angle of "bank 

angle of yaw 

a angle of attack: 

§ angle of sideslip (oov/Uq) 

6 angle -X axis makes with horizontal 

6 control setting, with appropriate subscript 
(6^ indicates that, values are in degrees 
and not in radians) 

X, 7, Z components of force along X., 7, and Z axes, 
respect ively • 

L rolling moment about X axis 

M pitching moment about Y axis 

IT yawing moment about Z axis 

E hinge moment 

b wing span 

c chord (of wing, unless otherwise subscripted) 
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S area (of wing, unless otherwise suhsc r ipt ed) 

I tail length (distance from center of gravity 
to tail post) 

m mass 

(jl relative density factor — S. 

l/2pSl 

inkv moment of inertia alsout X axis 



mkv^ moment of inertia ahotit Y axis 



mk^*" moment of inertia ahout Z axis 

t time 

D d/dt 

p density 

g gravity (32.2 ft/sec^) 

C^^ ro lling- moment coefficient ^ Ji--- — 

Vl/2pV^Sh/ 

Cm pi tching-moment coefficient / 

Vi/2pV"c/ 

Cj^ yawing-moment coefficient / ' II ^\ 

Cji hinge-moment coefficient / -? \ 

l,l/2pV^Sc; 

• Gy* liat eral-f orce coefficient f ^■ 

^ Al/2pV^S / 

Cj^ lift coefficient 

tail efficiency factor 

P empirical slipstream factor to account for con- 
tribution of fuselage, wing, and propeller to 
damping in pitch 
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lat eral-^ stability derivatives in terms of unit 
mass or moment of inertia of airplane. For 
example, 

^6 = 



P 
^6 

y 

^ longit-adinal-' stability derivatives in terms of 
unit mass or moment of inertia of airplane, 

X 

w For example, 

Xq 



Zw ^ m 



Zq 
Z6 



y 2, = ^-Z-Z35_ 



Za M., = MZM 



M 



mk 



Y 
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Subscripts are defined: 



solutions obtained when 61§ 



1 and 6Hg = 0 



U solutions obtained when 6L5 = 0 and SNg = 1 

a aileron 

e elevator 

r rudder 

t time 

0 initial value 

Primed symbols refer. to horizontal tail. 



SOLUTIONS OF THE LATERAL E<iUATIONS OP I'lOTIOM WHEN 
THE ROLLING DISTURBANCE IS UNITY AND TEE YAWIN& 
DISTURBANCE IS ZERO (5L5 = 1; 6H5 = O) 



The stability equation I*(D) =0 is obtained from 
equation (2) of the text by expanding the third-order de- 
terminant formed by the coefficients of the variables.. 
Thus, 



D-Lp 



-Lp 



F(D) = 







D-Np 
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III - le 



-L, 



-L 



-He 



= -/ + + iD^ 



U- (^r + 1-5;+ I^VLp LrNp + 



Let X^, Xg, X3, and "be the four roots of thii 

equation P(D) = 0 and form .the fo llowing p roduct s : 























- \) - 






(x^ -.X,) (x^ - 


^2) (^4"-* ^3 



) 



The solutions for any of the variables can be ex- 
pressed in terms of the products of the roots just formed 
and the stability derivatives. By use of the auxiliary 
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relat-ions tt- =.p and rr = r -so lut ions" -can be written 
a.z ■■ ■■ at 

directly for both 0;_,and' aI/. As it is Msually desirable 
in this work to detercine' ^. and if rather than p and 
r, it is convenient to solve for ^. ..and ^ direct Ij- and 
then to differentiate each of these solutions, respective- 
ly, to obtain p and r if these variables are wanted. 
The solutions follow: 



Angle of Bank 



?v. t 



A t 



A t 



A t 



T ® 



whe re 



A,- - A 



:U.-, 



9f 



9^ 



L_ = 



P. 



R 



7 ■ -^'rYfi 
+ Hp + ---C. 



S 



-■•^4T^r + 
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Angle of Sideslip Pj^ 

K^t A^t X^t 



where 



P, = J- -It 



'^^i Uo ft 



Uo H 



•^^3 Uo S 



ft. = J- (^4_:_Hr)_:_;_-£^ 



4 "0 T 



Angle of Yaw 



Xjt Xgt. At A^t 



where 
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S 

fj- 



L4 



vj/, -- 
Yawing Velocity r, 

^ dt 

The yawing velocity is most easily obtained "by direct 
different iat ion. 

Rolling Velocity 

The rolling velocity is most easily obtained 'by di- 
rect differentiation. 
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SOLUTIONS OP THE LATEEAl ECiUATIONS 05* MOTION WHEN 
THE YAWING DISTURBANCE IS UNITY AND THE EOLLING 
DISTURBANCE IS ZERO (8N5 = 1; 6Lg = O) 



Angle of Bank ^ 



N 



A,t A„t Agt X^t 



where 



■ ■ • „ ^N 

'0 "X "2 "3 4 



0 p 



= — 



frLR - Lr + >sLr 
= 
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irhere 



Angle of Sidesiip" p,^' 



X^t X^t X^t X^t 



No P 



gl>r 



- + 7\L^. - X. 



Uo IP 1 



Ug f 

R 



Angle of Taw .. 



where 
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"a 



3 



s 

-A 

1 






Uo " 


P Uo 














It. 




Uo 








S 










.3 P 


^§ 


g 


Uo 


Uo " 


^.^ u^ 






s 








It. 
Up 


^4^P 


Uo " 


u: 



Yawing Velocity r^j 
dt ^^N^ 

The yawing velbcity'is most easily obtained by direct 
dif f erent iat ion« • 



EolTing. Velocity 

The rolling velocity is most easily obtained by 
direct different iat ion. 
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Discussion OF TEE BCcUATION F(D) =0 

For most airplanes the solution, of FCD) = 0 will 
yield, both real and imaginary roots. In the solution of 
the equation, the real roo t s can first "be iso lat ed 
(Horner's method) and- extract ed f rom the equation. The 
imaginary roots can then be found by solving the result- 
ing quadratic. Because imaginary roots always occur in 
pairs, two of the roots will be of the form a '± ib, 
where a and b are constants* 

The components of the motion containing the imaginary 
roots can be combined conveniently into a single term in- 
volving only real numbers. ' In the case of the solu- 
tion is 



i i I Alt V. ^Vot' J X3t , X^t 

If and are the conjugate imaginary roots, 

^Jj ^ + e will also be imaginary. 

3 4 

Let X- = a + b; then will be imaginary and 

can be reduced to the form l^t-tZ. Further, 

. G + iH 

G + iE & -■ iH (jS ^ ga ... .. K. - .. 

and 

+ '= eat eos b (t + 7) 

where 
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Figure 8.- Variation of elevator angle required for trim with indicated airepeed. 
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